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ABSTRACT 

We present synthetic OH Zeeman splitting measurements of a super-Alfvenic molecular cloud model. We 
select dense cores from synthetic 13 CO maps computed from the largest simulation to date of supersonic and 
super-Alfvenic turbulence. The synthetic Zeeman splitting measurements in the cores yield a relation between 
the magnetic field strength, B, and the column density, N, in good agreement with the observations. The large 
scatter in B at a fixed value of N is partly due to intrinsic variations in the magnetic field strength from core to 
core. We also compute the relative mass-to-flux ratio between the center of the cores and their envelopes, 7Z^, 
and show that super-Alfvenic turbulence produces a significant scatter also in TZ^, including negative values 
(field reversal between core center and envelope). We find TZ^ < 1 for 70% of the cores, and TZ^ < for 12%. 
Of the cores with |Z?los| > 10 /iG, 81% have 1Z^ < 1. These predictions of the super-Alfvenic model are in 
stark contrast to the ambipolar drift model of core formation, where only 1Z^ > 1 is allowed. 

Subject headings: ISM: magnetic fields — stars:fromation — MHD — radiative transfer 
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1. INTRODUCTION 

The process of star formation in molecular clouds involves 
a complex interaction of turbulent velocity fields, gravita- 
tional forces, and magnetic fields. While the role of turbu- 
lence has been investigated in depth only in recent years, the 
importance of the magnetic field in star formation was rec- 
ognized long ago. On small scales, the magnetic field may 
support sub-critical prestellar core s against gravitation al col- 
lapse for an ambipolar drift time dLizano & Shulfl989h . and 
pro vide a way to shed angular momentum by magnetic break- 
ing {Mouschovias 1977, 1979) and by helping to power winds 
and jets (Pudrit z & Normanll 19831) . On larger scales, a strong 
enough magnetic field may even prevent the collapse of a gi- 
ant molecular cloud (GMC), pote ntially explain ing the low 
star-formation rate in the Galaxy (IShu et ai]|1987l) . 

As an alte rnative to this scenario of large scale mag- 
netic support, iPadoan & Nordlundl (1 19991) proposed a super- 
Alfvenic model of molecular clouds. In this model, the tur- 
bulent flows can compress the gas in all directions through 
shocks. Compressions with a component perpendicular to the 
magnetic field can locally increase the field strength, and the 
net result is a correlation of magnetic field and gas density, 
though with a large scatter. While in this model the mean 
magnetic field of a GMC is very low, its value inside dense 
cores and filaments (where it is usually observed) can be very 
large. Based on a number of tests comparing simulations with 
observations (IPadoan & Nordlundll 19991: IPadoan et al.ll2004l) . 
the super-Alfvenic model seems to explain the observations 
better than models assuming a strong magnetic field. 

The most direct way to infer the magnetic field strength in 
GMCs is by Zeeman splitting measurements. These measure- 
ments are very difficult, and are available in GMCs only for 
relatively dense cores, primarily from OH em ission lines (e.g. 
ICrutcher et al.ll993HTroland & Crutcherl2008l) . Only a hand- 

1 Helsinki University Observatory, P.O. Box 14, Tahtitorninmaki, FI- 
00014, University of Helsinki, Finland. 

- Department of Physics, University of California, San Diego, La Jolla, 
CA 92093-0424; ppadoan@ucsd.edu. 

3 Astronomical Observatory/Niels Bohr Institute, Juliane Maries Vej 30, 
DK-2100, Copenhagen, Denmark. 



ful of these measurements have provided detections of the 
field strength. Zeeman splitting is only sensitive to the line-of- 
sight component of the magnetic field averaged along the line 
of sight, and weighted by the OH emission (roughly speak- 
ing gas density and OH abundance). Thus, Zeeman split- 
ting does not directly provide the large-scale volume-averaged 
magnetic field strength, needed to understand the nature of the 
fragmentation process that leads to star formation. 

In this Letter, we relate the large-scale mean magnetic 
field to local Zeeman splitting measurements by a deductive 
method, predicting the observational properties of the theoret- 
ical model, rather than inducing the theory from the observa- 
tions. The super-Alfvenic molecular cloud model is obtained 
from a numerical experiment of three-dimensional, super- 
sonic, magneto-hydrodynamic (MHD) turbulence. Synthetic 
Zeeman splitting measurements are then simulated through 
radiative transfer calculations. This work is partly motivated 
by i) an improved interpretation of existing Zeeman splitting 
measurements by Crutcher (2008, pers. comm.), and ii) a new 
observational test of ambipolar drift by Crutc her et al.l (2008). 

A conclusion of the new analysis of Zeeman splitting mea- 
surements by Crutcher (2008, pers. comm.) is that the ob- 
served scatter in the relation between the absolute value of 
the line-of-sight magnetic field strength, |Blos|> an d the col- 
umn density, N, is partly due to a large scatter in the intrin- 
sic ma gnetic field strength, B, at any given N. As shown 
before (Pa doan & Nordlundl IT999). and as illustrated in the 
present work, this is a fundamental prediction of the super- 
Alfvenic model. If the mean field were very strong, its local 
value could not have a large scatter, and the observed scatter 
in |Z?los | could only arise as a result of random orientations, 
which seems to be ruled out by the new analysis of Crutcher 
(2008, pers. comm.). 

The new observ ational test of ambipolar drift by 
ICrutcher et al.l d2008l) . is based on the mass to flux ratio of 
the envelope relative to that of the core - henceforth re- 
ferred to as the "relative mass-to-flux ratio", 1Z^, as the am- 
bipolar drift model of core fo rmation predicts IZ^ > 1 (e.g. 
ICiolek & Mouschoviai (11994) -Fig. 3a). In this Letter we 
show that the super-Alfvenic model predicts a significant 
scatter in IZ^ (including negative values), with IZ^ < 1 in 
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FIG. 1. — Left: Simulated 13 CO (1-0) map of the model in the z-axis direction. The locations of the cloud cores are shown with squares. The circles indicate the 
locations of telescope beams used in the synthetic observations of three cores. Right: Line-of-sight magnetic field strength as calculated from Zeeman splitting. 



81% of the cores with |B L os| > 10 ^G. Values of 7\L M < 
1 in cores formed by turbulen t flows were first found by 
IVazqu ez-Semadeni et al. (2005), where a plausible mecha- 
nism for explaining that result was proposed. 

2. NUMERICAL SIMULATION OF SUPER-ALFVENIC TURBULENCE 

This work is based on a single snapshot of a supersonic 
and super- Alfvenic ideal MHD turbulence simula tion, run on 
a mes h of 1000 3 zones with the Stagger Code (Pado anet alJ 
120071) . We adopt periodic boundary conditions, isother- 
mal equation of state, random forcing in Fourier space at 
wavenumbers 1 <k<2(k=\ corresponds to the computa- 
tional box size), uniform initial density and magnetic field, 
random initial velocity field with power only at wavenumbers 
I <k <2. The rms sonic Mach number of the snapshot used 
in this work is A4 S = <7 v ,3d/c s = 8.91. 

The initial magnetic field is such that the initial value of the 
ratio of gas to magnetic pressure is /3; = 22.2. At the time 
corresponding to the snapshot used for this work, the rms 
magnetic field strength has been amplified by the turbulence, 
and the value of (3 defined with the rms magnetic pressure 
is (3 = 0.2. This corresponds to an rms Alfvenic Mach num- 
ber of .A/fa = (/3/2) 1 / 2 <Tv.3d/c s = 2.8, so the turbulence is still 
super- Alfvenic even with respect to the rms Alfven velocity. 
With respect to the Alfven velocity corresponding to the mean 
magnetic field = 22.2), the rms Alfvenic Mach number is 
much larger, M^i = (A/2) 1/2 ct v ,3d/c s = 29.7. 

3. RADIATIVE TRANSFER AND ZEEMAN SPLITTING 
MEASUREMENTS 

For the computation of synthetic Zeeman spectra the data 
cube is scaled to physical units. The length of the grid is 
fixed to L = 6 pc, the mean density to (n(H2)} = 150 cm -3 , 
and the kinetic temperature to T^ n = 20 K. The mean mag- 
netic field is (B X1 B,,B Z ) = (0.0,0.0,0.69) fiG and the rms field 
Brms = 6.45 /iG. We assume a constant fractional OH abun- 
dance of [OH]/[H] = 4.0 x 10" 8 recommended by ICrutcheri 
(fl979h . To compare the results with observations, the dis- 
tance to the cloud is fixed to 100 pc, yielding a total angular 
size of ^3.4°. The calculations are done for three orthogonal 
directions of a single snapshot. 

To simulate the selectio n of dense co res, we use our line ra- 
diative transfer program (lJuvelall 19971) to calculate 13 CO (1 - 



0) maps with an angular resolution of 50", after reducing the 
model resolution to 256 3 cells. Cores are handpicked in the 
13 CO maps for the Zeeman splitting analysis (the same anal- 
ysis on a set of cores selected with the Clumpfind algorithm 
(Willi ams et al.lll994l) yields practically the same results). To 
limit the contamination of envelope observations by emission 
from nearby structures unrelated to the core, we select only 
cores separated from each other by at least 10'. Some en- 
velopes can still be contaminated by weaker cores that were 
not selected for the analysis, but in calculations where we ig- 
nore the part of the envelope that is contaminated by weak 
cores we find no significant difference in the final results. We 
will discuss the core selection and the envelope contamination 
in detail in a forthcoming paper. 

The simulated OH Zeeman splitting observations are ob- 
ta ined with the same set of beams as used in the observations 
of ICrutcher et alJ d2008l) . The center of each core, i.e. the po- 
sition with the highest intensity in the CO map, is observed 
with a 3' (FWHM) beam, corresponding to measurements 
with the Arecibo telescope. Each of the envelopes of the cores 
is observed with four 8' beams, centered 6' north, south, east, 
and west of the core, simulating observations with the Green 
Bank Telescope. In the Zeeman splitting calculations we use 
a resolution corresponding to 512 3 computational cells. 

In our Zeeman splitting analysis we focus on the 
1665.40184 MHz OH line (hereafter the 1665 MHz OH line). 
The OH level populations are estimated assuming that the 
cloud is optically thin. The radiation field is approximately 
constant throughout the cloud, and the level populations de- 
pend only on the local density. Full radiative transfer calcula- 
tions with lower resolution models showed that the errors in 
the excitation temperature of the 1665 MHz OH transition are 
< 1 K. The coupled radiative transfer equations for the four 
Stokes parameters are integrated along the line of sight, tak- 
ing into account the effect of the magnetic field on emission 
and absorption, to obtain the observed I- and V-spectra. We 
use a bandwidth of 20 km s" 1 with 400 velocity channels. No 
noise is added to the synthetic spectra. 

We determine the line-of-sight magnetic field from the 
simulated observations by least-squares fitting the numerical 
derivative of the Stokes I-spectrum to the Stokes V-spectrum, 
as it is usually done with actual Zeeman splitting observations 
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FIG. 2. — i-e/?.' Estimated line-of-sight magnetic field strength as a function of hydrogen c olu mn density. The dots show the results from our simulated 
observations. The crosses show the results from observations listed in Troland & Crutcher 12008) and Crutcher 1 1999). Right: The distribution function for (N,B) 
from the simulations, using all the points in the three synthetic Zeeman splitting maps. The contours are at levels 0.01, 0.1, 0.4, and 1.5 (logjnCi^logjgCB)) , 



(e.g. ICrutcher et alJl!993t iBourke et al.ll2001h . Column den- 
sities of the cores are estimated from simulated OH spectra. 
Assuming that the line is optically thin, the column density of 
OH is obtained as 



jV(OH)[cnT 2 ] : 



4.04 x 10 1 



■W [Kkms -1 ] 



(1) 



where W is the integrated 1665 MHz OH line area, 7i, g is the 
background continuum brightness temperature ( 2.73 K), an d 
T ex is the excitation temperature of the transition (Goss 1968). 
We use T ex = 20 K in the analysis of our synthetic observa- 
tions. The left panel of Fig. 1 shows the simulated CO map 
and the locations of the selected cores in the z-axis direction. 
The right panel shows the corresponding line-of-sight mag- 
netic field as determined from the Zeeman splitting. 

4. RESULTS 

A total of 139 cores were selected from the 13 CO maps. 
The left panel of Fig 2 shows the absolute value of their line- 
of-sight magnetic field strength, |Blos|> as a function of their 
H2 column density, NQI2). Both the magnetic field and the 
column density are estimated from the simulated 1665 MHz 
OH observations with the 3' beam centered on the core. The 
figure also shows the results of 18 cm OH emission line Zee- 
man splitting observations from 45 cores that have been re- 
ported in the literature. T he observations include 34 cores 
from the recent survey of ITroland & Crutcher (2008]), and 
the 11 cores from the list in Crutcher (1993) that were ob- 
served in 18 cm OH emis sion lines, and were not included in 
ITroland & Crutcherl d2008l) . The observations and our simu- 
lations show a similar distribution of points in the (N, |Blos|) 
space. For a given column density, the inferred line-of-sight 
magnetic fields are widely scattered. This scatter in the super- 
Alfvenic model is partly due to intrinsic variations in the field 
strength from core to core, not only to the random orienta- 
tion of the field. The new analysis of Crutcher (2008, pers. 
comm.) finds that the distribution of observed |Blos| values 
implies a large intrinsic scatter of the field strength, as in the 
super- Alfvenic model. 

The upper envelope of the |Blos|-^V(H2) scatter plot is a 
useful constraint to compare the simulations with the obser- 
vations. However, given the small number of cores and their 
narrow range in column density, this upper envelope is not 
well defined. On the right panel of Fig. 2 we show the con- 
tour plot of |Blos| versus N(H2) for all the points in the 



three synthetic Zeeman splitting maps. In this case the up- 
per envelope is well defined, and approximately consistent 
with |fl LO s| ocJV(H 2 ) 2/3 . iPadoan & Nordlundl d 19991) found 
a similar upper envelope in the correlation between the three- 
dimensional distributions of magnetic field strength and gas 
density, B oc n 4 . This smaller exponent is consistent with 
the current result, because the gas density spans a much wider 
range of values than the column density. The power-law upper 
envelope derived from the right panel of Fig. 2 is replicated on 
the left panel, showing that it may be approximately valid also 
for the cores. 

Our simulation matches well the observed values of Blos 
in dense cores, despite the very low value of the mean field, 
(B) = 0.69 fiG parallel to the z axis, and (B) = 0.0 /iG in the 
other two orthogonal directions. If we compute the mean 
field by averaging over the Zeeman splitting maps, we ob- 
tain (Blos) =2.1 /iG in the z-axis direction (~ 0.1 /iG in the 
other two directions). Therefore, if one could ideally detect 
the Zeeman splitting everywhere on a cloud map like our syn- 
thetic map in the z-axis direction, the value of (Z?los)> aver- 
aged over the whole map, would be 2.1 /iG, while the cor- 
rect average value is only 0.69 /iG. More realistically, Zee- 
man splitting would be detected only in a few dense cores, 
and the mean value of those detections would overestimate 
the large-scale mean magnetic field by an even larger factor 
(by more than one order of magnitude in the case of the cores 
we have selected in the direction of the z axis). The volume 
average magnetic field is always overestimated by an average 
of Zeeman splitting detections because the measurements are 
essentially density weighted, and the magnetic field tends to 
be stronger at larger density, as a result of the B-n correla- 
tion responsible for the upper envelope of the Blos ~N scatter 
plot. Furthermore, only relatively large magnetic fields can be 
detected, biasing the average of Zeeman splitting detections 
to ward even larger valu e s . 

ICrutcher etalJ (2008) propose a new method for testing the 
ambipolar drift model of core formation. They compare the 
mass-to-flux ratio between the center of a core and its enve- 
lope, by measuring the relative mass-to-flux ratio, IZ^, de- 
fined as 



[N(H 2 )/B L0S ] 

C 



[jV(OH)/fi LOS ]c 



[jV(H 2 )/fi LO s]envelope LW(OH)/fi LOS ] envelope 



(2) 



The relative mass-to-flux ratio provides a direct test of 
the ambipolar drift model that requires TZ^ > 1 (e.g. 
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Fig. 3. — ie/?.' Relative mass-to-flux ratio for the selected cores as a function of column density. The red symbols indicate the cores with TZ^ < 0. The dots, 
crosses, triangles pointing down, triangles pointing up, and asterisks denote zero, one, two, three, or four field reversals in the envelope relative to the core center, 
respectively. Right: Relative mass-to-flux ratio as a function of inferred magnetic field strength in the central beam. The symbols have the same meaning as in 
the left panel. 

Ciolek & Mouschovias (1994) -Fig. 3a). This test avoids 
some of the problems faced by the direct measurement of the 
mass-to-flux ratio in the cores. For instance, it is not necessary 
to know the relative abundance of the tracer molecule, or the 
orientation of the magnetic field, although it is assumed that 
the relative abundance and field direction do not change from 
the core center to its envelope. We define the envelope A^OH) 
and Blos as the mean of their values from the four envelope 
beams. Because the Zeeman splitting gives also the sign of 
Blos, Hp is a signed quantity, with negative values indicating 
a reversal of the mean field in the envelope with respect to the 
central beam. Field reversals are interesting because they are 
not predicted by the ambipolar drift model. 

Figure 3 shows the absolute value of the relative mass-to- 
flux ratios from our simulation, plotted against the column 
density (left panel) and the absolute value of the line-of-sight 
magnetic field in the core center (right panel). The values 
of \TZfj,\ are widely scattered, with 91% of the cores in the 
range 0.2 < \1Z^\ < 5. While we do not find any clear trend 
with column density, there seems to be a correlation between 
|7^| and |fl L osl, with \1ZJ oc IB,™! -0 ' 26 . We find 1Z„ < 1 
for 70% 



'los |. wim |/t M | oc |fi LO s| 
of the cores, and TZ^ < for 12%. Of the cores 
with |Blos| > 10 /jG (values that could be observationally de- 
tected), 81% have IZ^ < 1. Even in cores with 1Z^ > 0, it 
is possible to find a field reversal in some of the four enve- 
lope beams, relative to the central one (i.e., the line-of-sight 
component of the magnetic field has opposing signs in the 
two beams). Figure 3 highlights these cores, with different 
symbols depending on the number of envelope beams with a 
reversed field direction. If we consider only the 36 cores with 
central values of |Blos| > 10 /iG, as weaker fields are less 
likely to be observationally detected, we find nine cores with 



one envelope field reversal, and one core with two reversals. 

5. CONCLUSIONS 

We have found that the super-Alfvenic model reproduces 
well the observed |Blos| ~N relation in dense cores, despite 
its very low mean magnetic field strength. Because the large 
scatter of |Blos| in the simulation is partly due to intrinsic 
variations of B from core to core, and not only to the effect of 
the random orientation, we interpret the observations to also 
be consistent with intrinsic variations in B from core to core. 
This interpretation of the observations, confirmed by the new 
analysis of Crutcher (2008, pers. comm.), suggests that the 
mean magnetic field on a larger scale cannot be as large as in 
the ambipolar drift model. 

The super-Alfvenic model also predicts that OH Zeeman 
splitting measurements of cores with central field values of 
|#los[ > 10 [iG should yield 1Z^ < 1 in more than 4/5 of the 
cores, and one or more field reversals in the envelope beams, 
relative to the central one, in 28% of the cores. In contrast, 
the ambipolar drift model of core formation predicts TZ^ > 1 
and no field reversal. If fu ture observational measurements of 
TZ^ dCrutcher et alj 2008) were to yield TZ p < 1 in most cores, 
as in the super-Alfvenic model, the ambipolar drift model of 
core formation would seem to be ruled out. 
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